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Abstract. We report on VLBA observations of a y-ray bright blazar NRAO 530 at multiple frequencies (5, 8, 15, 22, 39, 43, 
and 45 GHz) in 1997 and 1999. These multi-epoch multi-frequency high-resolution VLBI images exhibit a consistent core- 
' dominated morphology with a bending jet to the north of the core. The quasi-simultaneous data observed at five frequencies (5, 

^Nj ' 8, 15, 22 and 43 GHz) in February 1997 enable us to estimate the spectra of compact VLBI components in this highly variable 

source. Flat spectra are seen in central two components (A and B), and the most compact component A with the flattest spectral 
index at the south end is identified as the core. Based on the synchrotron cooling timescale argument, it is suggested that the 
^ I observed inverted spectrum of component C is caused by the free-free absorption (FFA), though the synchrotron self-absorption 

(SSA) model cannot not be definitely ruled out. While the SSA probably exists in component B, it is likely that the same FFA 
would produce the spectral turnover toward component B since the fitted FFA coefficients in both B and C components are 
almost the same. If so, the projected size of such an absorbing medium is at least about 25 pc. 
• By adding our new measurements to previous data, we obtain apparent velocities of two components (B and E) of 10.2 c 

' and 14.5 c, respectively. These are consistent with that the emergence of VLBI component is associated with the flux density 

, outburst, i.e. components B and E are related to strong y-ray flares in 1994.2-1994.6 and 1995.4-1995.5, respectively. We 

further investigate the spectral variability by making use of the single-dish measurements covering a complete outburst profile 
from mid-1994 to mid-1998. It shows a continuous increasing in the turnover frequency during the rising phase, and a gradual 
decreasing after passing the peak of the flare. Finally, we discuss the equipartition Doppler-factor (6eq) based on analysis of 
Q ' magnetic field and obtain S^gS of 3.7, 7.2 and 0.8 for components A, B and C, respectively, which are consistent with bf a larger 

flux density in component B, the non-detection of proper motion in component C and a bent jet. 

^ ' Key words, galaxies: quasars: individual: NRAO 530 - galaxies: jets - radio continuum: galaxies 
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^ 1. Introduction al. I1999> . To investigate the relationship between these flares 

d (at radio and y-ray) and the possible emergence of new compo- 
■ " ' NRAO 530 (PKS 1730-130) is an 18.5 mag QSO (Welch and ^^^^^^^ tj^i^ p^p^^. multi-epoch and multi-frequency 

SpinradT^) with a redshift z=0.902 (Junkkarinen[nHl. It vlbi study of NRAO 530. We will discuss the kinematics of 

IS variable at radio wavelengths (Medd et al.i972, Marscher the jet components, the possible absorption mechanism, and the 

and Marshall 1 1979> and is an emitter of X-ray (Marscher and spectral variability 

BrodericklWr). On the basis of its optical variability (Pollock Throughout this paper we adopt cosmological parameters 

et al. i_979) NRAO 530 is also classified as one of optically ^^^^7^ y.^ ^-i ^pc-i, Q^=0.27 and Qa=0.73. This results in 

violent variables (OVVs). ^ luminosity distance to NRAO 530 of 0^=5.8 Gpc, and an 

Figure □ (Top) shows the flux density history at 4.8, 8.0 angular-to-Unear size conversion factor of 7.8 pc mas"'. The 

and 14.5 GHz monitored by the Univers ity of M ichigan Radio spectral index a is defined by Sy oc v 
Astronomy Observatory (cf. Aller et al. ll985> . It can be seen 

that flux density changed greatly from 1994.5 to 1998.5. NRAO 3. VLBI Observations and Data Reduction 

530 underwent a dramatic millimeter flare beginning in 1994 

with two distinct flares nearly tripled the 90 GHz flux (Bower NRAO 530 was observed at 5, 8, 15, 22 and 43 GHz in 

et al. ll997> . EGRET reported increasing y-ray flux from the February 1997 with the VLBA plus a VLA antenna. It was 

direction of NRAO 530 between 1991 and 1995 (Hartman et also observed 3 times in 1999 at 39, 43 and 45 GHz with the 

VLBA. Detailed parameters of observations can be found in 



Send offprint requests to: S.-W. Feng, winfeng@shao.ac.cn Table [0 The data were recorded in VLBA format with a to- 

* e-mail: winf eng@shao .ac.cn tal bandwidth of 32 MHz (4 IFs) for left circular polarization 
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Fig. 1. Top: light curves of NRAO 530 at (A)4.8 GHz, (B)8.0 
GHz, (C)14.5 GHz from UMRAO (University of Michigan 
Radio Astronomy Observatory). Bottom: light curves of 
NRAO 530 observed during the period from 1994 to 1999 at 
(A)4.8 GHz, (B)8.0 GHz, (C)15.4 GHz from UMRAO and 
(D)22 GHz, (E)37 GHz from 13.7-meter antenna at Metsahovi 
Radio Observatory (Terasranta et al. '2004 l, arrows with letters 
a, b, c, d, e and f indicate different epochs plotted in Fig. 6. 



(LCP) (and also right circular polarization (RCP) in 1997). The 
observed data were cross-correlated with the VLBA correlator 
at Socorro, New Mexico of the USA. 

The post-correlation data reduction was carried out using 
the NRAO AIPS and Caltech DIFMAP packages. We first ap- 
pUed amplitude calibration (at frequencies >22 GHz atmo- 
spheric opacity corrections were made) using the gain curves 
and system temperatures measured at each station. Then fringe 



fitting was performed to determine the instrumental delays first, 
and further to solve for the residual delays, rates and phases. 
The data were then exported into DIFMAP for self-calibration 
imaging. These visibilities were time-averaged into 30-second 
bin with uncertainties estimated from the scatter of data points 
within the 30 s interval. The data were inspected for obviously 
bad points, most of which were near the beginning of scan 
when telescopes were still slewing to the sources. Images were 
made using the self-calibration and clean procedure for sev- 
eral iterations. A 1-Jy point-source model was employed at the 
start of each mapping process and the natural weighting of the 
data was used for all nine images shown in Fig.|2 The detailed 
paraments of images are summarized in Table|2] For a quantita- 
tive analysis of the structure in these images, we did model fits 
to the calibrated visibility data with elliptical Gaussian compo- 
nents. The results are listed in Table|3l A maximum uncertainty 
of 15% in the flux density was estimated from the uncertainties 
of the amplitude calibration and the formal errors of the model 
fits. The typical uncertainties of the component positions and 
their sizes were examined to be about 10% of the fitted val- 
ues. The brightness temperature (column 1 1 in Table|3j of each 
component was calculated using the following equation (i.e. 
Shen et al.[T997), 



Tb 



1.22X lO'^5v"^(0, 



ma J ^ mm 



r\i+z)K, 



(1) 



where S is the flux density in Jy, v the frequency in GHz, z the 
redshift, 6„,aj and fl,,,,,, sizes (FWHM) of the major and minor 
axes in mas, respectively. 



3. Results 

3.1. Components and Their Proper Motions 

The source shows different morphologies at different frequen- 
cies with more details in the core region revealed at higher fre- 
quencies (see Fig. |2j. The overall structure is dominated by a 
compact core plus a jet like feature extending approximately 5 
mas from the core along a position angle (P.A.) of 15° which 
further bends toward the north. Fig.|3lis the projected trajec- 
tory based on the locations of the VLBI components in NRAO 
530 observed in February 1997, inferring a continuously bend- 
ing path of the jet emission. Here, we assume the brightest and 
most compact component A (see Table 3) as the core because 
it is located at the south end of emission and has the flattest 
spectral index (see Sect. 3.2). The same core identification is 
also adopted by Jorstad et al. (2001 ). 

Component D at about 22-23 mas north of the core region 
was detected at 5 and 8 GHz (Fig.|2l(a) and (b)). This is proba- 
bly associated with the knot emission, at about 25 mas north of 
the core, seen from the 1.7 GHz VLBI observations (Bondi et 
al.^WE). At about 3-4 mas from the core with a P.A. of 10-15°, 
component C was consistently detected at frequencies of 5, 8, 
15 and 22 GHz (Fig. |2l(a), (b), (c) and (d)). This component 
was also seen by Tingay et al. <1998l l. Shen et al. iW^lt and 
Jorstad et al. i2001j . 

The 1997 simultaneous observations at 43, 22 and 15 GHz 
showed a consistent central core region (Fig.|2l(c), (d) and (e)) 
which can be fitted with 3 components. A, B and E (Table 3). 
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Table 1. VLBI Observations of NRAO 530 
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V 


epoch 


integration time 


bandwidth 


bit rate 


polarization(s) 


antennas 


(GHz) 


(year) 


(min) 


(MHz) 


(bit) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




1QQ7 10 


70 

Z.U 


32x2 


J 




VT RA4-VT A 1 




1 0Q7 1 






1 

1 


Rr'P T pp 


VT R A4-VT A 1 


15.4 


1997.12 


20 


32x2 


1 


RCP, LCP 


VLBA+VLAl 


22.2 


1997.12 


20 


32x2 


1 


RCP, LCP 


VLBA+VLAl 


43.2 


1997.12 


20 


32x2 


2 


RCP LCP 


VLBA+VLAl 


43.1 


1999.32 


17 


32 


2 


LCP 


VLBA 


43.1 


1999.40 


17 


32 


2 


LCP 


VLBA 


39.1 


1999.42 


19 


32 


2 


LCP 


VLBA 


45.1 


1999.42 


18 


32 


2 


LCP 


VLBA 



Notes: (1) Observing frequency in GHz; (2) Observing epoch in yr; (3) Total on-source time in minutes; (4) 
Recording bandwidth in MHz; (5) sampling rate in bit; (6) Polarization mode, and (7) Observing antennas. 

Table 2. Descriptions of VLBI images of NRAO 530 





V 


Epoch 


S peak 


Maj 


Min 


PA. 


Contours 




(GHz) 


(year) 


(Jy/beam) 


(mas) 


(mas) 


(deg) 


(mjy/beam) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(a) 


5.0 


1997.10 


5.7 


6.11 


2.21 


16.1 


2.9x(- 1 , 1 ,2,4,8, 16,32,64, 128,256) 


(b) 


8.4 


1997.10 


1.1 


3.80 


1.39 


15.5 


2.8x(- 1 , 1 ,2,4,8, 16,32,64, 128,256) 


(c) 


15.4 


1997.12 


8.0 


1.49 


0.49 


-8.0 


8.4x(- 1 , 1 ,2,4,8, 16,32,64, 128,256) 


(d) 


22.2 


1997.12 


6.8 


1.16 


0.32 


-9.9 


12.7x(-l, 1,2,4,8,16,32,64,128,256) 


(e) 


43.2 


1997.12 


4.5 


0.53 


0.19 


-7.3 


10.0x(-l, 1,2,4,8,16,32,64,128,256) 


(f) 


43.1 


1999.32 


1.7 


0.52 


0.19 


-8.2 


9.7x(-l, 1,2,4,8,16,32,64,128) 


(g) 


43.1 


1999.40 


1.8 


0.55 


0.19 


-9.0 


9.0x(-l, 1,2,4,8,16,32,64,128) 


(h) 


39.1 


1999.42 


2.0 


0.55 


0.20 


-7.0 


10.8x(-l, 1,2,4,8,16,32,64,128) 


(i) 


45.1 


1999.42 


1.7 


0.53 


0.19 


-8.3 


10.2x(- 1 , 1 ,2,4,8, 1 6,32,64, 128) 


Notes: (1) Image index (see Fig. 2); (2) Observing frequency; (3) Observing epoch; (4) Peak flux density; 



(6), (7) Parameters of the restoring Gaussian beam: the full width at half maximum (FWHM) of the major and 
minor axes and the position angle (P.A.) of the major axis. (8) Contour levels of the image. 



The fitted relative separations from A component of 0.35 and 
0.53 mas for components B and E, respectively at 43 GHz, 
are in a reasonable agreement with the results from 22 and 15 
GHz considering their >1 mas angular resolutions (See Table 
2). And it is also possible that the ~0.05-0.2 mas difference in 
the measured core separation for components B and E is due 
to the opacity shift, since the relative separation of VLBI com- 
ponent from the self-absorbed core is frequency dependent. At 
15 GHz, component E seemed too weak to be detected (also 
refer to Fig. |5] for its spectrum). These are comparable to the 
results obtained from four 43 GHz images at epochs 1995.79 
to 1996.90 (Jorstad et al. 2001 1 and 86 GHz image at epoch 
1995.38 (Bower et al. .1997i . Another four -43 GHz images 
(Fig. El (f), (g), (h) and (i)) made in April-May 1999 can be 
well fitted with 4 components: components A, B, E and a new 
component F at 0. 18 mas with RA. of 25°. 

By combining our new measurements of angular separa- 
tions with those published data in 1995 and 1996 (Jorstad et 
al.'200ri of components B and E, we obtain proper motions of 
0.21+0.02 and 0.30+0.02 mas yr"' (Fig.gJ, corresponding to 
superluminal velocities of 10.2+0.8 c and 14.5+1.1 c for com- 
ponents B and E, respectively. This infers the ejecting times for 
components B and E to be around 1 995 .46 and 1 994.76, respec- 



tively. These two epochs just fall in the outburst phase shown 
in Fig. 1 (also see discussion in Sect. 3.3). EGRET reported the 
increasing y-ray flux densities from the direction of NRAO 530 
during two periods (1994.2-1994.6 and 1995.4-1995.5) with 
the average fluxes of 64+26 pJy and 127+48 pJy (Hartman et 
aL 119991 . As EGRET could not observe NRAO 530 with dense 
time sampling and long integrations, the y-ray measurements 
have large uncertainties and can not rule out that y-ray flares 
were missed or that maximum of an ongoing flare was really 
seen. According to the existing data (VLBA and EGRET) of 
NRAO 530 we think that the emergence of components B and 
E are closely following the y-ray flares in NRAO 530. Activity 
of y-ray has been found to be correlated with infrared to X-ray 
activity in 3C 279 (Maraschi et al. 1994 1 and, to occur before or 
during a radio and millimeter outburst in many sources (Reich 
et al. .l993i) . 

For component C, we also plot its separation from the core 
(Fig-0 as a function of time measured in 1992, 1995, 1996 and 
1997 by Shen et al. HWTS, Jorstad et al. HM?), Bower and 
Backer (1998 1 and us (this work), respectively. The best fitted 
proper motion is about -0.01+0.13 mas yr corresponding 
to an apparent velocity of -0.35+6.04 c. This is quite differ- 
ent from the apparent velocity of 29+9/i 'c (Hubble constant 
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Ho=100 h km s ' Mpc deceleration parameter qo=0.1, and 
cosmological constant A = 0) reported by Jorstad et al. 42001 1 
(component D in their paper). Our new estimate with more time 
coverage in data indicates that there is no appreciable motion 
for component C. The previously reported large apparent ve- 
locity of 29+9/i^'c may be due to their limited epochs of data 
and thus the large uncertainty. This discrepancy may be related 
to the complex structure in component C (Fig. 2(c) and (d)), 
which results in a large error in the determination of its posi- 
tion. Such kind of diffused morphology could be due to internal 
jet structure, blending effects of shocks or some sort of inter- 
action. For the reason to be discussed in Sect. 4, we think it is 
mainly caused by the interaction of jet emission with the sur- 
rounding medium. 

3.2. Components' Spectra 

For the quasi-simultaneous observations in February 1997, the 
spectra of VLBI components (A, B, C and D) can be fit with 
Sy oc y " and the resultant optically thin spectral indexes 



are a^=0.21, aB=0.31+0.07, q'c=0.58+0.01 and q'd=0.91. 
Component A is regarded as the core of NRAO 530 since it has 
the flattest spectral index and it is at the south end of the over- 
all emission. Component E seems to show an inverted spectrum 
with a high turnover frequency (greater than 43 GHz) in 1997. 
This is likely to be due to the strong variability at 43 GHz as 
1999 data had a three times lower flux density. The fitted two- 
point spectral index is aE--l-5 between 22 and 43 GHz. 

In order to discuss the possibly dominant mechanism for 
the observed low frequency absorption in components B and 
C, both the synchrotron self-absorption (SSA) and the free-free 
absorption (FFA) are tried to fit the convex spectra with the 
following formulae, 

5v = 5ov"[l-exp(-T,v-"-2.5)] (2) 
for the SSA and, 

Sy = 5ov-"exp(-r/v-2') (3) 

for the FFA, where v is the observing frequency in GHz, Sq is 
the intrinsic flux density at 1 GHz in Jy, t/ and r, are the SSA 
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Table 3. NRAO 530 model parameters 
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V 


Epocti 


Component 


S 


r 





^ma j 


n 


P.A. 


Ti 


(GHz) 


(year) 




(Jy) 


(mas) 


(deg) 


(mas) 


(mas) 


(deg) 


(lO'^K) 


(1) (2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 



(a) 


5.0 


1997.10 


B 


5.60 








0.78 


0.23 


24.3 


2.9 








C 


0.60 


4.02 


14.7 


3.13 


1.70 


7.8 


l.OE-2 








D 


0.53 


21.92 


0.4 


16.30 


6.12 


13.2 


4.9E-4 


(b) 


8.4 


1997.10 


B 


7.73 








0.35 


0.21 


9.0 


3.5 








C 


0.64 


3.22 


12.3 


4.07 


1.65 


20.4 


3.1E-3 








D 


0.33 


22.89 


1.1 


13.60 


5.42 


13.9 


1.5E-4 


(c) 


15.4 


1997.12 


A 


0.91 








0.21 


0.21 




2.0E-1 








B 


8.01 


0.38 


16.5 


0.22 


0.10 


48.6 


3.66 








C 


0.46 


3.42 


10.2 


4.43 


1.70 


22.4 


5.9E-4 


(d) 


22.2 


1997.12 


A 


1.87 








0.26 


0.10 


12.0 


3.4E-1 








B 


6.73 


0.29 


24.5 


0.21 


0.16 


12.3 


9.4E-1 








E 


0.12 


0.73 


8.7 


0.14 


0.14 




2.9E-2 








C 


0.37 


3.24 


14.5 


4.89 


1.68 


13.1 


2.1E-4 


(e) 


43.2 


1997.12 


A 


1.62 








0.13 


0.04 


4.7 


3.9E-1 








B 


5.76 


0.35 


21.3 


0.17 


0.16 


-42.3 


2.6E-1 








E 


0.58 


0.53 


2.4 


0.66 


0.28 


4.3 


3.9E-3 


(f) 


43.1 


1999.32 


A 


1.63 








0.08 


0.05 


-0.3 


5.1E-1 








F 


0.34 


0.18 


13.4 


0.19 


0.19 




1.2E-2 








B 


0.28 


0.84 


27.5 


0.59 


0.13 


27.3 


4.6E-3 








E 


0.18 


1.42 


2.8 


1.11 


0.47 


17.7 


4.3E-4 


(g) 


43.1 


1999.40 


A 


1.86 








0.19 


0.04 


-12.4 


3.1E-1 








p 


0.22 


0.22 


20.8 


0.13 


0.13 




1.6E-2 








B 


0.26 


0.85 


27.9 


0.65 


0.18 


18.3 


2.8E-3 








E 


0.22 


1.33 


2.6 


1.22 


0.68 


19.5 


3.3E-4 


(h) 


39.1 


1999.42 


A 


1.78 








0.07 


0.04 


10.4 


9.6E-1 








F 


0.56 


0.17 


12.2 


0.34 


0.09 


45.5 


2.8E-2 








B 


0.31 


0.92 


26.0 


0.63 


0.20 


32.8 


3.7E-3 








E 


0.21 


1.48 


2.8 


1.24 


0.54 


23.6 


4.8E-4 


(i) 


45.1 


1999.42 


A 


1.65 








0.12 


0.02 


-11.2 


7.8E-1 








F 


0.16 


0.19 


18.5 


0.04 


0.04 




l.lE-1 








B 


0.34 


0.71 


27.4 


1.21 


0.12 


18.3 


2.7E-3 








E 


0.16 


1.29 


-1.6 


1.33 


0.21 


11.0 


6.5E-4 



Note: (1) Index of image (as shown in Fig.|2j; (2) Observing frequency; (3) Observing epoch; (4) Component 
in each image; (5), (6), (7), (8), (9), (10) Model parameters of Gaussian component: S=flux density, r=angular 
separation and (^=position angle of component with respect to the origin (0,0), full widths at half maximum 
(FWHM) of the major {9„aj) and minor (6„;„) axes and the position angle (P.A.) of the major axis; (11) 
Brightness temperature. 



and FFA coefficients at 1 GHz, respectively. We also performed 
the fits to the combined component (A+B) emission, which is 
believed to be represented by component B at low frequencies 
(5 and 8 GHz) due to the limited resolutions (see Table 3) and 
is resolved out at frequencies >15 GHz. All the fitting results 
are listed in Tablel^with the fitting curves shown in Fig.|5l 

Both the SSA and FFA fit the spectra of components B and 
C equally well. Thus, the spectral fit itself is not adequate to 
tell which mechanism is mainly responsible for the absorption 
towards components B and C. However, a short synchrotron 
cooling timescale seems to suggest that FFA plays an impor- 
tant role in component C (see Sect. 4). If so, that the FFA co- 
efficients for components B and C are almost the same (see 
Table 4) infers the presence of the FFA by the same suiTound- 
ing medium in component B. Then, one could estimate that the 



absorbing material should cover at least about 25 pc in size 
to produce the same FFA towards components B and C which 
were separated by about 3 mas in the sky in 1997. The interac- 
tion of the surrounding interstellar medium causing FFA with 
the jet emission in NRAO 530 can be used to explain the poorly 
determined proper motion in component C (see Sect. 3.1). 

For comparison, we also plotted the spectrum of the sum 
of all the VLBI components (labelled by T' in Fig.|5jl as well 
as the spectrum from single dish data (labelled by T in Fig.|5jl. 
The missing of flux density is persistently seen in all five VLBI 
frequencies. This is likely due to the effect of the resolution 
of the extended components. Therefore our amplitude calibra- 
tions, including the atmospheric opacity coiTections, should be 
reasonable. 
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Fig. 3. The projected trajectory of jet components observed in 1997: (Left) within central 4 mas, (Right) within 20 mas, showing 
a continuously bending jet moving away from the central core. A, B, E, C and D represent different VLBI components and the 
soUd Unes are connections of individual jet components. 



Table 4. Spectral fitting results with the SSA and FFA 



Parameters SSA FFA 





B 


C 


A+B 


B 


C 


A+B 


a 


0.31 


0.64 


0.14 


0.35 


0.79 


0.19 


SoOy) 


0.12 


0.02 


0.13 


21.77 


4.38 


15.59 


T 


148.31 


159.90 


101.40 


23.33 


20.88 


21.89 



3.3. Spectral Variability 

A complete outburst profile can be traced from mid- 1994 to 
mid-1998 (between two vertical lines in Fig.n(Top)). To an- 
alyze spectral variability during this radio outburst we select 
six epochs: 1994.42 (a), 1995.18 (b), 1996.39 (c), 1997.12 (d), 
1997.75 (e) and 1998.22 (f). Fig.glshows the spectral distribu- 
tion at these six epochs based on the data points in Fig. 1 (at 
frequencies < 37.0 GHz) and in Table 5 (at frequencies > 43 
GHz). It is quite clear in Fig.|6lthat the turnover frequency is 
increasing during the rising phase (represented by epochs (a), 
(b) and (c)) when the flux density is increasing (see Fig. 1). 



From 1997.12 (d) the turnover frequency starts to decrease 
when the total flux density decreases. This is because that high 
frequency photons lose energy quickly with jet expanding and 
thus the flux density at high frequency decrease too. In 1998.22 
(f) the spectrum returns to the level of quiescent phase. 

4. Discussion 

It is well known that the observed transverse velocity (JSapp - 
V apple, Vapp the apparent velocity) is related to the true ve- 
locity ifi) and the angle to the line of sight (9) by Papp - 
P sin 0/(1 -yS cos 0). Therefore, the maximum angle {ff) of 1 1 .2° 
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1992 1994 1996 1998 

Epoch (year) 



Fig. 4. Separation of components B, E and C from the assumed 
core component A as the function of the observing epoch. 

Table 5. Total Flux density data of NRAO 530 at v > 43 GHz 



V 


epoch 


flux 


error 


references 


(GHz) 


(Year) 


(Jy) 


(Jy) 




43 


1996.39 


16.6 


0.5 


Falcke et al.[T998l 


86 


1995.18 


11.2 


0.3 


Krichbaum et al.fT997l 


88 


1994.42 


7.8 


0.5 


Reuter et al.TWf] 


95 


1996.82 


12.6 


0.7 


Falcke et al., 19981 


107 


1996.39 


13.2 


0.4 


Falcke et al.[l998l 


142 


1994.42 


6.0 


0.5 


Reuter et al.[T997l 


150 


1996.82 


12.0 


1.2 


Falcke et al.[T998l 


215 


1995.18 


6.2 


1.1 


Krichbaum et al.ll997l 


375 


1997.12 


3.5 


0.5 


Robson et al.l2001l 


375 


1997.26 


2.7 


0.3 


Robson et al.l2001l 


375 


1997.75 


1.2 


0.1 


Robson et al.l2001l 



and 7.9° for components B and E can be estimated from their 
corresponding apparent velocities of 10.2 c and 14.5 c. 

Doppler factor (6) is defined as 5 = F"' ( 1 cos 0)"' , where 
F = (1 - ffiy^^^ is the Lorentz factor By means of the above 
expression of Papp, we can derive 5 = F + (F^ - Pipp ~ 1)''^- 
Then the Doppler factor of component B is 5 = F„„„ - 10.2 
when the Lorentz factor is F„„„ = (1 + Plpp)^^^- In the same 
way we compute the Doppler factor of another superluminal 
jet component E, 6 - 14.5. 

For another scenario in which no intrinsic acceleration or 
deceleration happened in the jet motion, i.e., jets are ejected 
with a same velocity /3 (or a constant Lorentz factor F), the 
changes of apparent velocities would be mainly due to differ- 
ent viewing angles (variable 6). The observed fastest apparent 
speed of 14.5c (for component E) in NRAO 530 will limit the 
constant Lorentz factor to be at least of 14.5, and thus we use 
Fmin - 15. Combining with the apparent velocity measure- 
ments, we can obtain the viewing angles of 9.7° and 1.6° for 
component B with Papp = 10.2c, and 4.9° and 3.0° for compo- 



SSA 




Frequency (GHz) 



Fig. 5. Components spectral fit of NRAO 530 in February 1997 
(top: spectra are fit with the SSA, bottom: spectra are fit with 
theFFA 

nent E with /3app - 14.5c. The corresponding Doppler factors 
are 4.0 and 26.0 for component B, and 11.3 and 18.7 for com- 
ponent E. 

Estimates of the magnetic field in a radio source can be ob- 
tained in two distinct ways. Firstly, by assuming the equiparti- 
tion of energy between the particles and the magnetic field, one 
can have (e.g., Pacholczvk i 1 970i> 

B'l = [4.5(1 -I- k)f(a, vuV2)LR-^f^'^ , (4) 

where B'^i in G is the magnetic field when the total energy of 
magnetic field and particles has a minimum value (equiparti- 
tion between the energy of the particles and the energy of the 
magnetic field), k is the energy ratio between the heavy par- 
ticles and the electrons (we used k=100), /(a, vi,V2) is tabu- 
lated in Appendix 2 of Pacholczyk (I1970> . a is the optically 
thin spectral index, vj and V2 are the lower and upper cutoff 
frequencies of the synchrotron spectrum with typical values of 
10^ and lO" Hz, R - is the radius of component in cm 
(here, is the angular diameter and d^, is luminosity distance), 
L is the synchrotron luminosity in erg s"' and can be expressed 
as L- ' Sydv (here, Sy is flux density at a frequency 

v). It should be mentioned that B'^'' is insensitive to the values 
of k and the fraction of the source's volume occupied by the 
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Fig. 6. Spectral variability seen in a flare from 1994.5 to 
1998.5. Flux densities at 43, 86, 88, 95, 107, 142, 150, 215 and 
375 GHz are listed in (Table|5}- Letters a, b, c, d, e, f represent 
different epochs indicated in Fig.[n 

magnetic field and the particles which is auusmed to be 1 . By 

XV-> 
" Sydv - S,„v„„ we can re-write 

B*"* ^4.0x 10"'*[/(a', 10^10")(1 -Hz)^c/^'0-^5,„v„,]2/7 ^ (5) 

where v,„ is the turnover frequency in GHz, S „, is the extrapo- 
lated optically thin flux density at the turnover frequency in Jy, 
and luminosity distance di and angular diameter are in Gpc 
and mas, respectively. 

Secondly, assuming that the inverted spectrum is due to 
SSA, one can estimate the magnetic field (i.e., Marscher ll983> . 



B^y" ^ lO-=b(a)&''^^„S;/5/(l + z) 



(6) 



where B^^" is the magnetic field in G, and b(a) is a tabulated pa- 
rameter dependent on the spectral index a (Table 1 in Marscher 
(Hlpl). For component B, q'=0.31 gives b(0.31)*3.0. 

From the SSA spectral fitting results for component B 
(0=V^ 



x6min-0-2 mas, Sm=8 Jy, v,„=10 GHz) and 6 =1 
(without consideration of relativistically beaming effect), we 



ma J 



obtained B^** 



580 mG and B, 



4 X 10^2 mG. The time 



scale for an electron to lose its energy through synchrotron 



■ 1 GHz 



This 



radiation is fi/2 - 2.76 x 10^ (yr) (f-c)" ( 
gives a cooling time scale of 1.1 x 10^ yr for component B, 
while the intrinsic cooling time scale in the source rest frame 
should be greater when Doppler factor is larger than (l-Hz). In 
the same way, we computed the magnetic field for component 
A, B'^ ^560 mG and Bj' 8 mG with S„,=1.87 Jy, v,„=22.2 
GHz, 0=0.16 mas, q'=0.21 and 6-1 since it seems to have a 
(presumably synchrotron self-absorbed) inverted spectrum at 
frequencies between 15 and 22 GHz (see Fig. 5). 

The ratio of the particle energy density (up) to the energy 
density of the magnetic field (m^) can be expressed as (i.e., 
Readhead rT994l ^ = i^^)"'"^- We get ^ ^ 6.9 x 10^ and 
— - 5.5 X 10'^ for component A and B, implying that both 
components A and B are particle dominated or that their emis- 
sions are relativistically beamed. In the former case, for com- 
ponent B this high ratio may account for the observed fast de- 
crease in the radio flux density (from 5.76 Jy to ~0.3 Jy) in two 



years even though its synchrotron cooling time could be much 
longer (^ 1.1x10^ yr). This is because that the high ratio of par- 
ticle to magnetic field energies may cause the inverse Compton 
catastrophe (Readhead ll994,) and, inverse Compton scattering 
of X-ray photons may further result in the y-ray emission (e.g., 
Hartmanet al. 1992). 

In the latter case, as the equations for the magnetic field 
(equations (5) and (6)) have different dependencies on the 
Doppler boosting factor, thus we can obtain 



6+i 



(7) 



where B*"* and B"'" are the measured magnetic field, B*^^, and 
B^l™ are the magnetic field estimated in the source rest frame. 
When B^'y ^ Bt'^", the Doppler factor in equation (7) is defined 
as the equipartition Doppler-factor, 6'''^. To derive equation (7), 
we have used the following relations, v°f = {j:;-;)v'm, - 
(f^)^^"5,'„ and 0°* = (■^)0' (Singal and Gopaf-Krishna ?), 
here0"*, v°f and5,^f are angular diameter, inverted frequency 
and corresponding ffux density in the observer's frame, respec- 
tively, and 0', v'„ and 5'„ in the source rest frame. Then we got 
the equipartition Doppler-factors for components A and B of 
S""^ ^3.7 and 6''^ - 7.2, respectively. These are comparable to 
the equipartition Doppler factors of 7.0 at 10 GHz (Marscher 
et al. 19771 and 5.2 at 15 GHz (Giiijosa and Dalv lf996t . 

Similarly, for component C we estimate B^* ^ 36 mG, 
B™ 29 X 10^^ mG, Uplu^ ^ 7.4 x 10"'^ 5'^ ^ 0.8 and 
the synchrotron cooling time of ~0.8 day, with the following 
parameters: S,„=0.65 Jy, v,„=7 GHz, 0=2.5 mas, q'=0.53 and 
^= 1 .0. As the magnetic field is strongly dependent on measured 
component parameters, the high magnetic field of component 
C may be due to its large size (~ 2.5 mas). Assuming that the 
equipartition is maintained in component C, we can estimate a 
size scale of 0.64 mas for component C so that both the cal- 
culated B"^ and B^"' are equal to 120 mG. The resultant syn- 
chrotron cooling time would be 7.6 yr, still too short for the 
SSA model. So, it is very likely that the FFA of the surround- 
ing medium dominates the observed absorption in component 
though we cannot completely rule out the SSA effect. 

As the Doppler-factors of jets are expected to increase with 
decreasing viewing angle for the assumed constant Lorentz fac- 
tor, the equipartition Doppler-factors also should decrease sys- 
tematically along a bent jet given that the viewing angle is in- 
creasing monotonically along the bending trajectory. This is 
probably the reason that the Seq of component C is smaller 
than the values of components A and B. However, the fact 
of S"^ > S'^ is consistent with that component B is brighter 
than component A in 1997. For component B, the Doppler- 
factor estimated from fiapp of 5 - T,„in - 10.2 is comparable 
to the equipartition Doppler-factor 6"^ - 7.2, suggesting that 
most likely jets have the minimum Lorentz factors in an equiv- 
alent state. For component C, there is no appreciable proper 
motion (Section 3.1), which suggests a Doppler-factor of ~1.0 

{5 - r„„„ = ^1 + 0lpp). This is consistent with the estimated 
equipartition Doppler-factor of 0.8, inferring that it approaches 
energy equivalent state between particle and magnetic field in 
the source rest frame. The fact that ^ 0.8 is smaller than 
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l+z= 1.902 for component C in NRAO 530 implies that it is 
Hobble flow that results in the observed departure from the en- 
ergy equipartition in component C. 

5. Summary 

We presented the high-resolution VLBI images of NRAO 530 
made at 5, 8, 15, 22, 39, 43, 45 GHz in 1997 and 1999. 
Combining our new measurements with the past data we esti- 
mated superluminal motions of components B and E at appar- 
ent velocities of 10.2 c and 14.5 c, respectively. These are con- 
sistent with that components B and E were ejected in 1995.46 
and 1994.76 following the y-ray flares of 127 and 64 pjy, re- 
spectively. Flat spectra are seen in central two components (A 
and B), and the most compact component A with the flattest 
spectral index at the south end is identified as the core. We de- 
rive the equipartition Doppler-factor based on the equations of 
(5) and (6), and obtain S'^'^s of 3.7, 7.2 and 0.8 for components 
A, B and C, respectively. These are consistent with a larger 
flux density in component B, the non-detection of proper mo- 
tion in component C and a bent jet. The synchrotron cooling 
timescale argues that the FFA is responsible for the inverted 
spectrum of component C though the SSA efl'ect can not be 
completely ruled out yet. While for component B both the SSA 
and FFA could play a role, we propose that the FFA may be 
at work given that almost the same FFA coeflicients were fit- 
ted for both components B and C. If so, the absorbing material 
would have a projected size of at least 25 pc. 
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